Edge effects in the two-dimensional spin-| Heisenberg antiferromagnet 
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We use quantum Monte Carlo simulations to study effects of free edges in the two-dimensional 
spin- 1 Heisenberg antiferromagnet. We find that the magnetic response of an edge is smaller than 
the bulk susceptibility. This counter-intuitive quantum effect can be traced to enhanced antiferro- 
magnetic nearest-neighbor spin correlations, i.e., tendency to local singlet formation, at and close 
to the edge. These correlations form a comb-like pattern, which can be reproduced with a simple 
variational valence-bond state. We also study rough edges, and find that these instead significantly 
enhance the susceptibility, due to local sublattice imbalance impeding singlet formation. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Cx, 75.40.Mg 



Among the intricacies of strongly-correlated quantum 
systems, the roles of various defects, such as impurities 
and boundaries, are intriguing. Defects inevitably ef- 
fect experiments, to an extent often not precisely known. 
On the other hand, they can also serve as useful exper- 
imental probes of correlated quantum states [H, 0, H IH • 
Theoretically, impurity effects in quantum antiferromag- 
nets have been studied extensively in one @, 0] and two 
dimensions @, @, H, EH EH, EH- Impurities can effec- 
tively cut isolated spin chains into finite segments with 
free ends, which leads to particularly strong deviations 
from bulk magnetic properties in one dimension [5|. A 
good example is the quasi-one-dimensional antiferromag- 
net Sr2Cu03, for which the NMR line exhibits a broad 
background [2fl which is well accounted for by the lo- 
cal magnetic susceptibility (Knight-shift) distribution of 
open-end Heisenberg chains @, IE]- In two dimensions, 
free edges can be expected to have less dramatic conse- 
quences, because of the typically small ratio of bound- 
ary to bulk, and not much attention have been paid to 
them. However, with the increasing focus on nano-scale 
materials, the boundary physics should become accessi- 
ble (or unavoidable, depending on the perspective) also 
in two-dimensional antiferromagnets. It is therefore im- 
portant to establish what edge effects to expect based 
on prototypical model systems such as the Heisenberg 
hamiltonian. In this Letter we take some steps in this 
direction. 

We use the approximation-free stochastic series expan- 
sion (SSE) quantum Monte Carlo method [l3| to study 
the antiferromagnetic (J > 0) Heisenberg hamiltonian, 
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(1) 



where Si are the usual S = | spin operators and 
denotes nearest-neighbors on a square L x L lattice. We 
consider systems with completely open boundaries as well 
as semi-open ones, which are periodic in one direction 
and open in the other direction. The absence of corners 
and the translational symmetry along the open edges in 



the semi-open systems allow easier access to an infinite 
edge. On the other hand, it is also interesting to study 
corner effects in the fully open systems. Experimentally, 
it is likely that typical samples would have some rough- 
ness, leading to effects not captured by the smooth edges 
of the L x L lattices. We therefore also study systems 
with irregular edges, constructed according to a scheme 
described further below. We are interested in the mag- 
netic response of the edges, and will also investigate how 
this is related to changes in the spin-spin correlations 
relative to those in fully periodic systems. 

In analogy with the impurity susceptibility previously 
considered for systems with isolated vacancies or added 
spins [^, E3| we define an edge susceptibility, 
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Xa ~ X0 
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where Xa is the total magnetic susceptibility for a sys- 
tem with a free edges; a = for periodic systems, a = 2 
for semi-open boundaries, and a = 4 for fully open sys- 
tems. In all cases there are N = L 2 spins and the total 
susceptibilities are given by 
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Xa = j;(M 2 z ), M z 
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The normalization with aL in ^ reflects the natural 
assumption that the difference in response should scale 
with the total length of the edges of the open or semi- 
open systems. One would intuitively anticipate xe > 0, 
as the edge spins should be free to fluctuate more than 
those in the bulk. Surprisingly, this actually does not 
hold when T« J and L — ► oo. Fig. [1] shows the tem- 
perature dependence of xe for semi-open systems. The 
scaling of xi — Xo with the edge length is confirmed as 
the xe curves for different L coincide for increasing L. 
For T ^S> J, the spins contribute independently (4T) _1 
to the susceptibilities and, consequently, xe vanishes. In 
the limit T/J — > 0, we must have xe —> for any (even) 
L because the ground state is a singlet regardless of the 
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FIG. 1: (Color online) Edge susceptibility of semi-open sys- 
tems. Error bars are smaller than the symbols. The line is a 
log-lin fit to the size-converged data. 

boundary condition (i.e., M z = 0), which is seen explic- 
itly for L = 4, and 8. Focusing on the L — > oo converged 
data, decreasing T initially leads to an increasing xe, in 
line with the expectation of X2 > Xo due to enhanced 
fluctuations of the edge spins. However, a maximum is 
reached at T/J 0.5, below which xe decreases and 
becomes negative. The temperature dependence below 
T/J w 0.1 is consistent with a logarithmic divergence; 
Xe oc — ln(J/T). This behavior can be contrasted with 
the single-impurity susceptibility, which, for both vacan- 
cies and added spins, is always positive when T — > 0, 
diverging as a Curie form with a log correction [^, |T(| • 
In the case of the edge problem considered here, where 
the defect is non-magnetic (the number of spins is the 
same in the periodic and open systems), there is no a 
priori reason to expect a divergent xe- We have also run 
simulations for the classical Heisenberg antiferromagnet 
with semi-open boundaries and in that case find xe to 
converge to a positive constant as T — > 0. Hence, the neg- 
ative divergent edge susceptibility has to be attributed to 
quantum effects beyond those included in the "renormal- 
ized classical" description 1J] of the Heisenberg model. 



When normalized by the total number of spins L , in- 
stead of the edge length L, the divergent edge suscepti- 
bility would be a negligible correction to the bulk sus- 
ceptibility, which is constant as T — + 14 1 . However, in 
principle a local susceptibility can be accessed in NMR 
experiments through the Knight shift 0, [Hj], and, pro- 
vided that sufficient sensitivity can be achieved and the 
edges are smooth enough (both of which may clearly pose 
challenges), the divergence should be detectable. We thus 
also study the spatially resolved susceptibility, defined 
for a site % as Xa(i) = 0{SfM z ). For a periodic system, 
a = 0, there is no dependence on the location i, whereas 
in the semi-open system X2(*) depends only on the dis- 
tance R of i from the edge; in either case J2i Xa{i) = Xa- 
We define the position dependent edge susceptibility, 



Xe(R)=X2(R)-Xo/L 2 , 



(4) 



which obeys J2rXe{R) — Xe- Fig- El shows this quan- 
tity at two different temperatures at which the results are 




FIG. 2: (Color online) Position-dependent edge susceptibility 
at distance R from the free edge of a 64 x 64 semi-open lattice 
at two different temperatures. The temperature convergence 
for R = 0, 1 is shown in the inset. 

size-converged for L = 64. Here it is seen that a large con- 
tribution to the negative edge susceptibility comes from 
the second line of spins from the edge, R = 1. The re- 
sponse of the edge line, R = 0, is always larger than 
the bulk value, however, and xe(0) seems to vanish, or 
become very small, as T — > 0. Beyond R w 10, xe(R) 
becomes very small, indistinguishable from zero, also for 
the lower temperature. The available data suggest that 
Xe(R) should become negative for any R at sufficiently 
low T (and L — ► oo). Considering the log-divergent Xe, 
we should have xe(-R) oc when T — > 0. 

We will next show that the negative edge susceptibil- 
ity is related to enhanced local spin correlations close 
to the edge. The spin correlators (Si • Sj) for nearest- 
neighbor sites i,j (bonds) form a pattern of weak and 
strong bonds. These correlations provide a measure of 
the amplitude of spins i,j forming a singlet. In addition 
to calculating the correlations with the SSE method, we 
have also used a variational state in the valence-bond 
basis, from which some additional insights are gained. 

A valance-bond basis state for TV spins is a product of 
N/2 singlets (a, b) = (talb - laU)/V^, where a and b 
are sites on different sublattices of the bipartite square 
lattice. Any singlet state \^>) can be expanded in this 
over-complete basis; 



l*)=E^)l( a l' 5 l)---(^/2'^/2)) I 



(5) 



where v € {1, . . . ^ !} labels the different bond configu- 
rations. In the amplitude-product state of Liang et al. 
[ill], the wave function coefficients are products of real 
amplitudes h(a,b); 



N/2 

i>{v) = n He 



(6) 



For a periodic system the amplitudes depend only on the 
bond lengths (the x and y separations of the two sites) , 
h(a^,b^) = h(Xi,yl), but in an open system they depend 
on both site coordinates (up to reflection and rotation 
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FIG. 3: Bond patterns at an edge of an L = 16 system ob- 
tained with (a) an optimized amplitude-product state (almost 
exact), (b) a state with amplitudes h(r) = e~ r , and (c) clas- 
sical dimers. The line widths correspond to —(Si • Sj} in the 
range [0.315,0.451] and [0.316,0.466] in (a) and (b), respec- 
tively, and average dimer occupation € [0.167,0.500] in (c). 



symmetries). It is known that the periodic Heisenberg 
model can be very well described by this simple state 
3 12 1- The amplitudes decay as h(r) 
is the bond length 
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r , where r 
To study boundary (including 
corner) effects, we have optimized all amplitudes for a 
16 x 16 fully open lattice, using the optimization method 



discussed in 17| . Also in this case the amplitude-product 
state provides a very good description of the system, 
with the energy deviating by less than 0.1% from the 
approximation- free result obtained using the SSE method 
at very low T. The bond pattern is also almost identi- 
cal in the variational and SSE calculations; the result is 
shown in Fig. GDJa). The correlations are the strongest 
at the corners, and at the edge they form a comb-like 
pattern. This comb-structure is repeated on alternating 
columns away from the boundary, with a rapidly decay- 
ing amplitude, as shown in Fig. |4] based on SSE calcula- 
tions. The strongest bonds are enhanced by more than 
10% compared to the bulk. This enhancement is associ- 
ated with higher amplitudes for local singlets within the 
combs, which is clearly consistent with the reduced edge 
susceptibility. The correlation modulations are seen in 
Fig. |4]to decay more rapidly with R than the edge suscep- 
tibility in Fig. [21 On the other hand, the nearest-neighbor 
correlations represent only the dominant contribution to 
local singlet formation, and thus the two properties are 
not easily related to each other quantitatively. 

In order to elucidate the origin of the comb-like bond 
pattern, as well as the plaquette structure at the cor- 
ners, we have calculated the correlations also in a non- 
optimal state where the bond-amplitudes are very short- 
ranged, h(r) = e~ r , which corresponds to a spin liquid 
[l6j |. In spite of this state being very different from the 
actual Neel ground state of the Heisenberg model, a very 
similar bond pattern forms at the corners and edges of 
the open lattice, as shown in Fig. EJb). This indicates 
that the edge pattern is rather insensitive to the long- 
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FIG. 4: (Color online) Deviation AC tj = C, 
nearest-neighbor correlation Cy = | (Si ■ Sj 
value Cfj versus the distance R from an edge of an L — 64 
semi-open system at T — J/32. Integer and half-integer R 
correspond to bonds ij parallel and perpendicular to the edge, 
respectively. The inset shows magnified R > 2 data. 

range correlations of the state — it is essentially governed 
by the hard-core nature of short valence bonds. To fur- 
ther illustrate this point, we show in Fig. [3fc) the bond- 
occupation pattern of the classical dimer model (with 
only short bonds, averaging over bond configurations us- 
ing Monte Carlo sampling) . Here there is more of a ten- 
dency to plaquette formation at the edge, which, how- 
ever, changes into a uniform comb-structure away from 
the corners of larger lattices. These results show that the 
gross features of the boundary correlation pattern of the 
Heisenberg model is dominated by the entropy of short 
valence bonds, and thus we argue that this is at the heart 
also of the reduced edge susceptibility. 

The site dependent susceptibility Xa(i) can be used 
to derive experimental consequences, such as the NMR 
line-shape 
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However, it is unclear whether the 
edge effects we have discussed so far could be observed 
experimentally. Samples consisting of extremely small 
fragments are most likely required to distinguish any edge 
features from the NMR bulk signal, and the fragments 
appearing in powders hardly have long smooth edges; 
more likely they have irregular shapes. In principle one 
could carry out simulations for a suitable ensemble of 
clusters. However, not knowing the actual structure of 
clusters that could be expected in experiments, we here 
consider a simple model for roughness added to the Lx L 
open systems discussed above, with the aim of studying 
the robustness of the smooth-edge effects. 

Our roughness model amounts to traversing the 4(L — 
1) boundary sites of an L x L lattice and removing each 
spin with probability p or coupling a new spin to it with 
probability p (and doing nothing with probability 1 — 2p) . 
In order to have the same number of spins in the peri- 
odic Lx L systems and these rough-boundary system, we 
only study samples with the same number of added and 
removed spins, and, furthermore, we only consider clus- 
ters with equal numbers of spins on both sublattices (so 
that the total ground-state spin remains 0) . The impurity 
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FIG. 5: Edge susceptibility for systems with rough edges. The 
line is a fit to the size-converged intermediate- T data. 

susceptibility xe, Eq- ©, then still vanishes for finite L 
both in the limits T — > and T — > oo. Figure [5] shows re- 
sults averaged over several hundred random boundaries 
for each L. Here p = 1/3, corresponding to maximum 
roughness. We again observe a non-trivial logarithmic 
divergence, but, in contrast with the smooth edge, xe 
is always positive. The prefactor is about 5 times larger 
than in the smooth-edge case. Thus we conclude that the 
roughness has completely changed the nature of the edge 
effect. For less rough boundaries (smaller p) the prefac- 
tor of the log divergence is reduced, and for some very 
small p we expect to recover the negative factor pertain- 
ing for the smooth boundaries. We have not yet carried 
out systematic studies of this, however. 

Within the picture of the reduced susceptibility of the 
smooth edge arising from local singlets, the different re- 
sponse of the rough edges can be understood as an im- 
peded singlet (short valence bond) formation due to lo- 
cal sublattice imbalance. This has been previously ex- 
amined in diluted systems (randomly placed vacancies) 
18]. In that the percolation point is approached 

local sublattce imbalance leads to localized magnetic mo- 
ments, which interact and form low-lying states at an 
energy scale below the normally lowest-lying "quantum 
rotor" states of the Heisenberg antiferromagnet on finite 
clusters. Here, for the rough-edge problem, we have only 
studied the static magnetic response, but it would clearly 
be interesting to study also other aspects of the rough 
boundaries, e.g., their excitations. 

In summary, we have found that smooth open edges 
in the two-dimensional S = \ Heisenberg model have a 
smaller magnetic response than the bulk, contrary to the 
naively expected enhancement of fluctuations due to the 
smaller number of neighbors of the edge spins. We have 
explained this surprising effect in terms of local singlet 
formation at the edges, which in turn can be regarded as 
a consequence of entropy maximization in a valence-bond 
description of the system. In sharp contrast to smooth 
edges, rough boundaries lead to an enhanced magnetic 
susceptibility. We have argued that this is due to local 
sublattice imbalance ("dangling spins"), which impedes 
local singlet (short valance bond) formation. For both 



smooth and rough boundaries, the edge susceptibility of 
an infinite system diverges logarithmically as T — > 0. 
These studies also demonstrate that edges should have 
profound effects on the magnetic response of nano-scale 
clusters, and that details of the boundary texture are 
important. 

The smooth-edge effects that we have pointed out here 
have very recently been examined using field-theoretical 
methods by Metlitski and Sachdev [19|- The negative 
edge susceptibility originates from low-lying spin waves. 
The prefactor of the log-divergence, the slope in Fig. [H is 
in reasonable agreement with the prediction of [191 ] . The 
comb-structure in the edge correlations was argued to 
be a short-distance phenomenon, as we have also shown 
here, beyond the standard 0(3) continuum filed theory 
description. It can be understood in terms of proximity 
to a phase transition into a valence-bond-solid state. 
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